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Abstract 

G\ ' 

£\\ • Refraction and diffraction of waves in natural crystals and artificial crystals formed by 

anisotropically scattering centers are considered. A detailed study of the electromagnetic wave 
c/2 ' refraction in a two-dimensional photonic crystal formed by parallel threads is given by way of 

example. The expression is derived for the effective amplitude of wave scattering by a thread 
{~i !■ (in a crystal) for the case when scattering by a single thread in a vacuum is anisotropic. It is 

established that for a wave with orthogonal polarization, unlike a wave with parallel polariza- 
tion, the index of refraction in crystals built from metallic threads can be greater than unity, 
• , and Vavilov-Chrernkov radiation becomes possible in them. The set of equations describing the 

J>^| dynamical diffraction of waves in crystals is derived for the case when scattering by a single 

■ center in a vacuum is anisotropic. 
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Because a most general approach is applied to the description of the scattering process, 
the results thus obtained are valid for a wide range of cases without being restricted to either 
electromagnetic waves or crystals built from threads. 



2 Introduction 

^ I Creation of metamaterials has recently become an area of vigorous research worldwide. The so-called 
CN ; electromagnetic (photonic) crystals built from, e.g., metallic split rings or parallel metallic threads [I], 
including threads with dimensions within the nanometer range [2] are actively being studied. Such 
photonic crystals can be used, and are already being used, for solving various tasks, particularly 
^ ! in antenna microwave technology [Tj. In addition, crystals built from periodically strained parallel 
^ ! metallic threads can serve as resonators in volume free electron lasers ( VFEL) (3J HJ E] . 

The interaction of electromagnetic waves with photonic crystals is accompanied by the phenomena 
of refraction and diffraction. As is known, the refractive index of the medium formed by randomly 
distributed scatterers is related to the amplitude of scattering by a single center as follows [5] 



n 2 = l + -j£A{0), (1) 
where p is the density of scatterers, A(0) is the amplitude of forward scattering. 
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However, it is shown in [3 [8] that if scatterers are located periodically (e.g. crystal), then for a 
correct description of the refraction process, one should use a modified expression for n 2 
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where Q3 is the volume of the unit cell of the crystal. If scattering by a single centers is elastic, 
formula fl2]), in contrast to ([I]), leads to a physically correct result: the imaginary part of the refractive 
index equals zero. Note that fl2]) is derived under the assumption that scattering by a single center is 
isotropic, i.e., the amplitude A is independent of the scattering angle. In the present paper, refraction 
of waves in crystals built from anisotropically scattering centers is considered. Metallic threads are 
the simplest example of such scatterers: scattering by such threads of electromagnetic waves with 
the electric-field vector polarized orthogonally to the thread axis is anisotropic for all wavelengths 

PUD]. 

This paper is arranged as follows: Section 1 gives a detailed analysis of a nonplane waves scattering 
by a single thread. Section 2 describes the method used for finding the refractive index by the 
example of crystals formed by isotropic scatterers. Section 3 considers the case of crystals formed by 
anisotropic scatterers. 



1 Scattering of Electromagnetic Waves with Orthogonal Po- 
larization by Thread 

1.1 Plane Wave 

Before we start to consider scattering of a cylindrical wave by a thread (cylinder), let us recall the 
case of a plane wave. Assume, as usual, that the radius R of the thread is much smaller than its 
length L, and so in analytical treatment, the thread can be considered infinitely long. 

The solution to the problem of diffraction of a plane electromagnetic wave by an infinite cylinder 
can be found in the form of a series over Bessel and Hankel functions [11] . Let a wave with perpen- 
dicular polarization (vector E is perpendicular to the axis of a cylinder) be scattered by a cylinder 
placed in a vacuum (FigJT]). The axis of the cylinder coincides with the z-axis of the rectangular 
coordinate system. Let us also introduce a cylindrical coordinate system (r, <p, z), as is shown in 
Figffl 




Figure 1: Diffraction of a plane electromagnetic wave by a cylinder. 
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For simplicity, we shall consider the case when the wave vector k is perpendicular to the axis of 
the cylinder. Then the field of the scattered wave can be written in the form 



He = e z £ i n ciH n (h 

(3) 



n=— oo 



oo 

—in 



4 = - E 



n=— oo 



e r H n (kr) - e^kH'^kr) 



r 



e -inip 



where the coefficients are calculated by the formula 

± -J n {k 2 R)J' n (kR) + ±J' n {k 2 R)J n {kR) 



J n (k 2 R)Hi iy (kR) - ±J' n {k 2 R)H£\kR) 



(4) 



Here J n is the Bessel function of order n, H n is the Hankel function of the first kind of order n, 
k 2 = ky/e, e is the dielectric permittivity of the thread (for metals, e = 1 + Anai/uj, where u is 
the frequency of the electromagnetic field, a is the conductivity, and the magnetic permittivity \i is 
assumed to be equal to 1), vectors e r , e v , e z are the unit vectors of the cylindrical coordinate system, 
and the amplitude of the incident wave is assumed to be equal to 1. 

Considering the limits of Bessel and Hankel functions of small argument, one can find that for 
a perfectly conducting cylinder, the relationships Cq « — and ^> hold for all |n| > 1 at 
kR < 1. If a cylinder is not a perfect conductor, then the first equality is violated, but the absolute 
values of the coefficients q} and c± x remain comparable. What is more, if the dielectric permittivity 
of a cylinder is e — 1 ~ 1, then q} c^. Thus, at kR <C 1, in it suffices to take account of only 
the terms n = and n — ±1. Note here that in considering a wave with parallel polarization in the 
long-wave limit, cj} 3> c| for all n^O, and so in the series for the field, one can take account of the 
term n = alone. 

Returning to our problem, let is recall that the two equations ([3]) are related through Maxwell's 
equations, and so the second equation can readily be obtained from the first one by formula E e = 
— TotH e . With thus eliminated electric field and with due account of the above remarks concerning 
the values of the coefficients c n (it is further assumed that kR 1), one can write the scattered 
wave in the form: 

H e = e z ^ sc = e z {c H (kr) + 2ic x H x {kr) cosy?} , (5) 

where the superscript _L on the expansion coefficient is dropped. 

If a plane wave H Q = ^$e z = e lkx e z , is scattered by a cylinder placed at the origin of coordinates, 
then the total wave field is expresses as a sum of the field of the incident and scattered waves, i.e., 

$ = ^ + \£ sc = e ikx + c H (kr) + 2ic l H x {kr) cos (p. (6) 

Using the asymptotic expression for Hankel functions of large argument and the integral repre- 
sentation of Hankel functions, the above expression can be presented in the following form, provided 
kr > 1: 

00 ■, r 

nr-^ ^ (7) 

— oo 

where A{ip) = — ^(c + 2c x cosy?) = A + A x cos (p. 

Similarly to a three-dimensional case, by A(p) one should understand the amplitude of scattering 
of an electromagnetic wave by a thread at an angle (p [12] . 
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It should be noted that in contrast to a diverging spherical wave that characterizes scattering in 
a three-dimensional case, in a two-dimensional case, a diverging cylindrical wave is formed. 

In view of (jH]), one can readily write the expression for the field in the case when the axis of a 
cylinder lies not in the origin of coordinates, but at point r\ 



= # + = e ikr + e im ■ A inH (k\r- f x |) - e ik?1 ■ A^H^r - r x \) cos(fc,r - r t ). (8) 



1.2 Cylindrical wave 

Let now a cylindrical wave \l/o = Ho(kr), diverging from the origin of coordinates, be incident 
onto this scatterer, which is placed at point f±. To analyze the scattering process in this case, one 
can decompose a cylindrical wave into elementary plane waves, which are scattered according to a 
well-known law (181) . For the Hankel function, in particular, one can use the following representation: 

oo oo 
I r gikr J r & ikr 

H (kr) = - / dk v = - / dk x . (9) 



-oo — oo 



For the purposes of this paper, we do not need to know the exact decomposition, suffice it to 
know that the wave \l/o = Ho(kr) can be presented as a sum of plane waves as follows: Ho(kr) = 
J f(k,ky)e lkr dk y , or H Q (kr) = J f(k, k x )e lkr dk x , where k = \k\ — const and f(k,k y ) is a certain 
function. 

Thus, by decomposing the initial wave and in view of (JBJ), one can immediately write the expression 
for a scattered wave 



®sc = J f{k,k y y k ^A Q mH Q {k\r-r 1 \)dk y - j f{k,k y )e lk ^A 1 TxH 1 {k\r-r 1 \)cos{k,r-r 1 )dk y . (10) 



-oo 



Because A 0jl and if ,i(^l^ — are on ly dependent on the absolute value of the wave vector k and 
independent of k y , they can be removed from the integration sign. Then, according to (jHJ), the first 
integral in fTTOj) equals Ho(kri), and f TTOj) can be rewritten in the form: 

= H (kn) ■ A nrH (k\r- fi|) - A^H^f- r x \) ■ h, (11) 

where 

oo 

J 1= ! f(k,k y )e ik " 1 cos(k,f-f 1 )dk y . (12) 



1 In (fT0|) . e l ^ Tl = e l V^^v Xl+lk y Vl ■ at fc y > k, this quantity is a damped plane wave whose wave number is k y , but 
not k. At large values of k y (k y R 3> 1), the coefficients c n (n > 1 ) for such a wave are comparable with Co and C\ 
(see (01). And so, in integrating over the domain of large values of k y , it would be desirable to add in ([TU)) the terms 
of the form 

H n (k y \r - r x \) x / c„(k y )cos(na)dk y , 

J ijkl - fc 2 

where a is the angle between vectors k and r — r\. Simple estimates, however, show that the absolute values of these 
integrals are small (~ e~ k y Xl <C 1), and so the corresponding additions can be neglected. 
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To evaluate the integral I\ , let us represent the cosine of the angle between vectors k and r — T\ in 

k(f—T\) 
k\r—f\ | 



the form: cos(fc, r — f\) = jf^-^l- Then we have 



oo 



h = ■ - / k y )e fkfl dk y = 

k\r — 7*1 1 7T 7 

— oo 

oo 

= - ^^~ r ' -grad! / /(A;, k y )^dk y = -i T ~ T \ • gra^ ff (*ri), (13) 

— 00 

where the subscript "1" on the gradient symbol means that differentiation is performed with respect 
to the coordinates of the point f\. The differentiation yields grad x ifo(^ r i) = ~Hi(kr 1 )^-, and then 

7*1(7** — 7*1) 

ii = iHi(kri) — — ; — — - = iHiikri) cosm, r — fi). (14) 
fi\r — 7*1 1 

The above yields that the total wave field for scattering of a cylindrical wave has the form: 

* = # (/cr) + A i7ri?o(A;ri) ■ #o(&|r - ri|) - Ai7rzifi(A;ri) • i?i(A*|r- fi|) cos(fi,r - Fi). (15) 

At large distances from the origin of coordinates (kr\ 3> 1), one can use the serial expansion iHi(z) ~ 
Hq(z) (l + t-). Substitution of this expression into (fl5]l gives 



* = HQ{kr)+AQmHQ{kr 1 )-H Q {k\r-r 1 \)-A X 'KH Q {kr 1 ) ^1 + -i?i(A;|r— ri]) cos(ri, r-ri). (16) 

Comparing this equation with (jSJ), one can conclude that with increasing distance between the origin 
of coordinates (the central point of the diverging cylindrical wave) and the scatterer, the difference 
between the amplitudes of scattering of cylindrical and plane waves by a thread diminishes, as might 
be expected. 

Let us briefly consider the case when the initial cylindrical wave has the form \l/ — Hi{kr) cos (p. 
By reasoning along the same lines, we obtain similar formulas, except for the expression for the 
integral I\. Now it equals 

I' 1 = -i- rT z, — zrr ■ gradi Hi(kri) cos^i- (17) 
k\r — 7*1 1 

Now differentiation yields the following expression 

{2 1 f 1 1 
kH (kri) H\{kri) > cos^i • 1 Hi(kri)e x , (18) 
7*1 J 7*1 7*1 

where e x is the unit vector of the x-axis. The expression for a scattered wave will now have the form: 

* = ^0 + ^sc = H x {kr) cos <£> + {AoH^kr^ cos^i} inH (k\r- f*i|)- 
Ai compx (-iH {kr 1 ) + ^-H^kr^J j 7ri?i(A;|r — ri|) cos(f 1 ,f - 7*1) — 

— H\{kr x ) \ -KH x {k\r- fi|) cos(f- ri,e x )- (19) 

KTi J 

A particular case when the scatterer is placed on the y-axis, i.e., <f% = ±|, seems to be of interest. 
Here we have cos<£>i = 0, and the expression for a scattered wave simplifies appreciably: 

^sc = i ^H 1 (kr 1 )H 1 (k\f- fi|)cos(f- f u e x ). (20) 



It should be noticed that here the scattered wave does not vanish despite the zero amplitude of 
the incident wave at the scatterer location: \&o(^l) = -Hi(fcri) cosy?i = 0. Though there is not any 
paradox because at <pi = ±|, only the magnetic-field vector equals zero, H = ^oe z , while the 
electric-field vector has the form: E = ^rotH = — ^Hi(kri)e y , i.e., is nonzero. Thus, when 
(pi = ±|, only the electrical component of the initial wave is scattered by the thread. 

2 Refraction of Waves in Photonic Crystals Formed by Isotrop- 
ically Scattering Elements 

Let us consider an infinite, two-dimensional crystal composed of periodically arranged scatterers. A 
well-known example of such a crystal is a photonic crystal built from parallel metallic threads [3], 
Fig. El 
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Figure 2: Crystal made of parallel metallic threads 

For definiteness (but without loss of generality!), we shall consider such a crystal. Let the coordi- 
nates of the threads in the lattice be (x mn , y mn ) = (ma, nb), where a and b are the lattice spacings, m 
and n are the integers; the axes of the treads are parallel to the z-axis of the rectangular coordinate 
system. 

— # 

One should distinguish between the cases when vector E of the incident wave is polarized parallel 
to the threads and when it is polarized perpendicular to them. In the first case, scattering at kR <C 1 
is isotropic and A(ip) = A Q , while in the second case, scattering is anisotropic even at kR <C 1, and 
the amplitude has the above described angular dependence: A = A + A\ cos tp. Certainly, when 
kR > 1, in expression ([6]) one needs to retain more terms in the expansion ([3]), and the angular 
dependence of the amplitude can take a more complicated form. 

In this section, we shall assume that scattering by the centers is isotropic with the amplitude A . 
Let us suppose that an electromagnetic wave ~ e iqx propagates in a crystal in the positive direction 
of the x-axis. @ This wave results from the summation of diverging cylindrical waves i%H (k\r — f mn \) 
radiated by all the threads. The amplitude of these waves has the form $ = ^ e iqXrn ", where $o 

2 The wave vectors k and q are assumed to be perpendicular to the z-axis. The case of arbitrary incidence of a wave 
onto the z-axis will be considered individually. 
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is independent of the position of the thread in the crystal. Using the method described in [T3] , let 
us find the relation between the wave numbers k and q. Let us assume that the local field at the 
location of each thread is a sum of all waves coming from all other threads. Particularly, the local 
filed acting on the thread placed at the origin of coordinates (m = n = 0) can be written as follows: 

* loc = $ iir Yl e UtXm »H (kr nhn ) = $ in Yl e^ m ^ (V(am) 2 + (M 2 ) • (21) 

(m,n)^(0,0) (m,n)^(0,0) 

Since scattering is isotropic, all these waves are scattered by this thread with a known ampli- 
tude, equal to the amplitude A of scattering of a plane wave by the thread, producing a diverging 
cylindrical wave with the amplitude $o : 

$ = a v 1oc = A m% Y eiqamH ° (V(H 2 + N 2 ) • 

(m,n)^(0,0) 

We thus come to the dispersion equation for finding q 

_1_ = e tqam H (ky/(am) 2 + (bn) 2 ) % = S(k,q). (22) 

U (m,n)^(0,0) 

The sum in (|22|) is taken as described in [T3], and here we only give the resulting expression: 

J 2 / kb \ 2 sinka 2^ 1 sinki^a 

S(k,q) = -l-t-<- log — + C }+— : + 7 > -7-r 

I 7T \ Ait I hn rns hn — nns art h ^— i l„\ n ) 



7T V 4vr y A;6 cos ka - cos ga 6 ^ ^ cos fci n) a - cos qa ^\n\ 

(23) 

where k^ = i^J (2irn/b) 2 — k 2 , C ~ 0.5772 is the Euler constant. Let us assume that the scattering 
amplitude is sufficiently small: |Ao(/c)| "C 1 and the refractive index is close to unity: \n — 1| = 
\q/k-l\ < 1. 

One can easily demonstrate that in this case, the second term between the braces is a dominating 
term, and far from the diffraction conditions (ka ^ irn) it can be presented in the form: 

2 sinka 4 1 4 1 . . 

(24) 



kb cos ka — cos qa k 2 ba n 2 — 1 k 2 Q.2 n 2 — l' 

where Q 2 = &b is the crystal unit cell area. We thus come to the following dispersion equation: 

1 4i 1 

-1 - 



Aoin k 2 VL2 n 2 — 1 ' 

from this we get 

9 47T An , . 

"= 1 + «5'TTlk- (25) 

At this point, it should be mentioned that the same result was obtained in [3] in a different way. 
In the case of elastic scattering (e.g. perfectly conducting threads), formula (j25j) gives a real value of 
the refractive index. This is easily verified by the optical theorem. It will be recalled that the optical 
theorem relates the imaginary part of the forward scattering amplitude to the total scattering cross 
section: ImA(O) = For a two-dimensional case, the differential scattering cross section is related 
to the scattering amplitude as follows: ^ = ^^(y?)! 2 . 

One can arrive at the same result using a slightly different method: start with the consideration 
of scattering of a plane wave by a one-dimensional grating built from threads and then proceed to 
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Figure 3: Diffraction of a plane wave by a one-dimensional array of scatterers 

the case of infinite crystals. As this method is more convenient for the purposes of the analysis given 
in the following section, let us also describe it here. 

Let a plane wave be scattered by a one- dimensional array of scatterers, as is shown in Fig. El 
The scattered wave in this case is a sum of cylindrical waves of the same amplitude F, which are 
radiated by each thread 



(x, y) = F inH (k\/x 2 + (y — bn) 2 ). 



But the amplitude F of the plane-wave scattering by a thread in the presence of other scatterers 
differs from the amplitude A Q of scattering by a single thread. With the help of the general methods 
used for describing multiple scattering [6], the amplitude F can be expressed in terms of Aq as follows: 



F = A + A Fj2^ H o(kb\n\). 



(26) 



The physical meaning of equations (126]) is obvious: there are two waves scattered by each thread 
(particularly, by the thread placed at the origin of coordinates): the plane wave e lkx and the wave 
= F ^ i7rHo(kb\n\) scattered by all other threads with the amplitude Aq. The sum fl26l) is taken 

in HT 



S 1 = Y,Mkb\ 



n\ 



2 

kb 



i- ( 

7T 



log 



kb 

Air 



2 1 



1 



y/(2irn/b) 2 -k 2 27r|n| 



(27) 



It is still assumed that Aq <§C 1, then using (I26|) and (j27[) . one can derive the following expressiorH 
for the effective scattering amplitude: 

A 



1 — inSiAQ 1 + zttAq 



; 2n a 



3 To simplify the derived expressions, in writing (|28|l and (|29[) . we have made the inessential assumption that kb < 2ir. 
Quite simple, but cumbersome calculations (subsequent consideration of the cases 2ir < kb < An, An < kb < 67r, and 
so on) show that this assumption has no influence on the final result, which remains valid for large values of k (far 
from the diffraction conditions) too. 
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Thus the wave field resulting from scattering of a plane wave by a one-dimensional grating has the 
form: 

00 2 ?7T 

■qr = e ikx + F inH {ky/x 2 + (y - fen) 2 ) « e fe + — Fe* |x| , (29) 

n=— 00 

where F is defined from (12 8p . and is a sum of the wave fields of the plane wave incident onto 
the grating and the wave scattered by the grating. The scattered wave propagates in the positive 
direction of the x-axis at x > 0, and in the negative direction at x < 0, its amplitude in both cases 
being ^rrF. This quantity can be treated as the "amplitude of scattering" of a plane wave by a 
grating built from threads. 

Now let us consider the propagation of waves in a crystal formed by an infinite periodic system of 
plane gratings regularly spaced at an interval a and composed of threads. Use the same method as 
before to derive the dispersion equation. The wave incident onto the grating placed at the origin of 
coordinates is a sum of waves scattered by all other gratings. Assuming that the amplitude of these 
waves has the form $ = Q e iqam with $ being independent of the coordinates due to the periodicity 
of crystals, one can write the following expression 

kb ^— ' 

from which one can derive the dispersion equation (compare with f|22l ) 

I — p e iqam e ika\m\ ^ 

kb 

The sum in (l30j) can easily be calculated using the geometric progression sum formula, and is equal 
to 

e igam e ika\m\ = _j _ • Sm ka _ _j _ ^_ . 1 

^— ' cos ka — cos qa ka n 2 — 1 

m^O ^ 

By substituting the obtained value into equation (I5UI) and using (EHjl . we get the already known 
result 

2 ! , 4tt 

Pfi 2 1 + inA ' 

Note here that similar consideration of a three-dimensional crystal case gives a well-known for- 
mula (El). 



3 Refraction of Waves in Crystals. The Case of Anisotropic 
Scatterers 

Assume now that the scattering amplitude has the form A = A + A\ cos p, i.e., scattering is 
anisotropic. Let us use the approach described in the previous section to consider the propaga- 
tion of plane waves in a crystal composed of anisotropic scatterers. 

Let a plane wave having the amplitude e lkx be incident onto a one-dimensional grating composed 
of anisotropically scattering centers (Fig. |3]). Knowing how a cylindrical wave is scattered by such 
scatterers (see (TTol) . (IT91) ). one can determine that the angular dependence of the amplitude F{ip) of 
scattering by each thread in the presence of other threads is the same as that of the amplitude A(ip) 
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of scattering by a single thread F(<p) = Fq + F\ cos tp. The system of equations for finding Fq and 
Fi has the form (see (1151) and (El): 



F = A Q + A Q F J2^ H o(kb\m\), (31) 
Fi = A 1 + A 1 F 1 ^mH x {kb\m\)-{kb\m\y 1 . (32) 

The sum appearing in equation f[3"2l can be calculated using the known value of the first sum f[2"Tl) 
and the following equality for the Hankel function: 



Hx{kb\n\) 2z 1 

fcfe|n| 7ik 2 b 2 n 2 k 2 

Integration gives: 



kH (kb\n\)dk. 



o 



_ v ^ ffi(fc6|m|) _ 2 l_i_( r ]_ kb 2tt 2 \ 2z v ^ / 27r|n| fc 2 fe / 47r 2 r2 2 

2 ~ V kb\m\ kb~2~n ( 2 + ° g 4^ + 3fc¥j"B^ I _ 6 4dd ~ V ~P 

(33) 

Then, as before (see fT2"8l) ). we can write the expressions for the amplitudes F and Fi using equations 
(J2ZD and (J33D 

F n = — 1° h (34) 

l-mS 2 1 + i\A\ - %%A\ ' 1 ; 

So the scattered wave field at a point with coordinates (x, y) has the form: 

V(x,y) = e ikx + A jr ittH (k\J (y - nbf + x 2 ) +A X ^ (-tt) -H^lkyjiy- nbf + x 2 ") cos p 

n=— oo ^ 7 n=— oo ^ 7 

(36) 

where 

cos^ n = ^== X (37) 

Summation of these series is given, e.g., in [9]; finally, for sufficiently large distances \x\ from the 
grating, we have 

■q f = e ^ + ^l(F ±F 1 )e ik{x{ , (38) 
kb 

where the sign "+" refers to the case when x > and the sign "— ", to the case when x < (forward 
and backward scattering, respectively). 

Now let us consider a crystal composed of a number of such gratings regularly spaced at an 
interval a. The wave which falls onto the grating placed at the origin of coordinates, is a sum 
of plane waves scattered by all other gratings. Let us assume that the amplitude of these waves is 
$ = ($o i $i)e*' ?am , where <3>o and $i are independent of the grating number due to the periodicity 
of the crystal. Then the wave has the form 

— — --$ e iqam e ika\m\ _| "Hl^ J ^.qam e ika\m\ _ gigam ^,ka\m\ I ^ ^39^ 

m^O k.m=— 00 »i = l J 



/ 1 ■ 
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In view of fl39|) and (1381) . one can write the following system of equations for the amplitudes 

*° = F °{lr* oS3 + ^r* iS4 )' (40) 
= F >{lr*« s < + lr*'4 (41) 

where the sums S3 and S4 are equal to 

1 

n 2 — 1 ' 
1 

n 2 - 1' 

By equating to zero the determinant of the system (|40|) - (|4ip . one obtains the dispersion equation of 
the form: 

1 = ^ tf> + Fx) S 3 + (^) 2 W (52 - S 3 2 ) . (42) 

Simple (though cumbersome) arithmetic transforms of this equation with substituted values of the 
sums S3 and S 4 as well as the expressions for scattering amplitudes F and F x (formulas (1341) - (155]) ) 
give the final expression for the refractive index of a crystal 

9 47T f An Ai 1 

"- 1+ i^'{ra + TTipr}- (43) 

Note that a three-dimensional case can be considered in a similar manner. As a result, if the 
amplitude of scattering by a single center has the form A{9) = A + Ai cos 9, then the refractive 
index of the crystal can be written as follows: 

2 Air j A At 1 

71 " 1+ F^\l + ^o + l + z|A 1 J ' 1 j 

where f2 3 is the volume of the crystal unit cell. 

Under the assumption of smallness of the scattering amplitude and with the above-selected form 
of its angular dependence, the obtained expressions hold for any values of k far from the diffrac- 
tion conditions. Moreover, the scatterers can obviously be arbitrary, not necessarily threads. In a 
particular case of crystals built from metallic threads, substitution into (j4*3|) of An and A\ for the 
wave with perpendicular polarization (see formulas (T4|), (J7D) gives for the refractive index the value 
ttR 2 

n 2 = 1 + — — > 1, which means that in such crystals the Vavilov-Cherenkov effect can be observed 

PEED]. 

In quantum mechanics, the scattering process is usually described using the transition matrix T 
(which is non-Hermitian in general). Recall that the diagonal element of this matrix is proportional 
to the amplitude of forward scattering. Equation ()44j) reflects the fact that for an ordered medium 
(crystal), in the expression for the refractive index, the diagonal element of the matrix T must be 
replaced by the diagonal element of the Hermitian reaction matrix K (the properties of the matrix 
K see in [To]). 



S 3 = e iqam e ika\m\ _ _ ;, Sm ka ^ _j _ 2l 

/ -— ' cos fca — cos oa fca 

sin ga 2i 



S 4 = y^ie 



-iqam iqam\ ikam 



cos fca — cos qa ka 
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A detailed analysis shows that in view of the above results, the equations describing the dynamical 
diffraction of waves in crystals must be modified. Their general form, of course, does not change: 



: o / 



(45) 



$(f) = ^ v {k + r)e l{U ^\ 



(46) 



T 



where g(f) is the structure amplitude, f is the reciprocal lattice vector of the crystal. However, in 
the expression for the structure amplitude, the amplitude of scattering by a single thread is replaced 
by the effective amplitude of wave scattering by a thread in the crystal: 



This result agrees well with that given in [3 [8] for the isotropic case. According to [8], when 
calculating the structure amplitude, one must exclude from the imaginary part of the scattering 
amplitude the contribution to the total cross section that comes from elastic coherent scattering. 

Conclusion 

This paper considers the process of propagation of waves in natural crystals and artificial crystals 
formed by anisotropically scattering centers. The interaction of waves with individual scatterers is 
described in terms of the scattering amplitude. Special consideration is given to taking account of 
multiple rescattering of the initial wave by the centers in a crystal. The obtained expression ( F43]h 
which relates the refractive index of a crystal to the scattering amplitude, differs from the known 
formula for the refractive index of non-periodic media ([T]) and allows one to correctly describe wave 
attenuation in crystals. In particular, if scattering by a single center is elastic, the refractive index of 
the crystal, calculated in accordance with (T4"5|) . is a real value, e.i., the attenuation is absent, whereas 
application of a conventional formula ([T]) in this case leads to an erroneous conclusion that n has a 
nonzero imaginary part. 

The equations derived in this paper, which describe the field in crystals, coincide with standard 
equations of the dynamical diffraction theory in crystals [161 C3 EH]- However, in the expression for 
the structure amplitude g(f), the amplitude of scattering by a single center must be replaced by 
the effective amplitude of wave scattering by a center located in the crystal (which is described by 
the Hermitian reaction matrix K). This enables one to correctly describe the effect related to the 
attenuation of coherent waves in crystals. 

In the present paper, a detailed consideration of the electromagnetic wave refraction in a two- 
dimensional photonic crystal built from parallel metallic thread is given by way of example. An 
interesting result for such a crystal is that the index of refraction for a wave with E polarized 
perpendicular to the threads is greater than 1, and the Vavilov-Cherenkov effect can be observed 
in the crystal [HI [10]. A general approach applied here to the description of scattering enables one 
to obtain the results that are valid for a wide range of cases without being restricted to either 
electromagnetic waves or crystals built from threads. They can be of interest, in particular, for 
studying diffraction of cold neutrons in crystals, investigating of various nano crystalline materials, 
designing metamaterials with prescribed properties, etc. 




(47) 
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